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We find an interference effect for electron-plionon interac- 
tions in coupled semiconductor quantum dots that can domi- 
nate the nonlinear transport properties even for temperatures 
close to zero. The intradot electron tunnel process leads to a 
'shake up' of the phonon system and is dominated by a double- 
slit-like interference effect of spontaneously emitted phonons. 
The effect is closely related to subradiance of photons in a 
laser-trapped two-ion system and explains the oscillations in 
the nonlinear current-voltage characteristics of coupled dots 
observed recently. 

PACS: 73.23.Hk, 71.38.-Fi, 42.50 Fx 



Spontaneous emission is one of the fundamental con- 
cepts of quantum mechanics that can be traced back to 
such early works as that of Albert Einstein jl]. An ex- 
cited state of a single atom decays exponentially due to 
the coupling to photons. In a system of two atoms inter- 
acting via the common photon field, the decay splits into 
a sub- and a superradiant channel. This effect is a pre- 
cursor of the famous Dicke superradiance phenomenon ||] 
and was verified experimentally in the spontaneous emis- 
sion of photons from two trapped ions only three years 
ago§. 

Recently, in a completely different physical system, the 
emission of phonons from two artificial atoms has been 
observed Q]. Here, the coupling to the phonon degrees of 
freedom turned out to dominate the non-linear electron 
transport through semiconductor double quantum dots 
even at mK temperatures. 

Double quantum dots are well-defined artificial sys- 
tems for the study of interaction |5j and coherent time- 
dependent Q effects. Here, we propose a theory for a 
new non-linear transport effect in double quantum dots 
which corresponds to the Dicke effect, i.e. the collective 
decay of real atoms. In our theory, the tunneling of single 
electrons through coupled artificial atoms is renormalized 
by the interaction with piezoelectric acoustic phonons 
and leads to an orthogonality catastrophe of the phonon 
bath if an electron tunnels between the dots. This 'boson 
shake up' effect is determined by an effective den- 
sity of states p{uj) of the phonon modes Q that couple to 
the tunnel process. These interfere like in a double slit 
experiment when interacting with the electron densities 
in the two dots. As a result, p{lo) shows oscillations on a 
scale LOd '■= Cs/d, where Cs is the speed of sound and d the 
distance between the centers of the two dots. It turns out 
that the non-linear current peak as a function of the dif- 



ference £ between the two relevant many-particle energies 
is determined by the shape of p(uj = e/K). Furthermore, 
this quantity is analogous to the rate for emission of sub- 
radiant photons from two laser-trapped ions |^ , when Cg 
is replaced by the speed of light and d by the distance of 
the ions. Thus, both phenomena are physically closely 
related. This provides the microscopic mechanism for 
the oscillations observed recently in a double dot current 
spectrum We predict that future experiments with 
artificial atoms can exploit this analogy to real atoms in 
more detail. In particular, coherent effects such as su- 
perradiance can be manipulated by gate-voltages and 
external leads in such systems ||lo| ]. 

In our model, we consider a double quantum dot con- 
sisting of a left and a right dot (L and R, respectively) 
which are connected via a tunnel barrier. Each dot is 
connected to an electron reservoir in thermal equilib- 
rium with chemical potentials pL (source) and pn < pL 
(drain). For the physical phenomena we are interested 
in, it is sufficient to restrict the dot Hilbert space to 
the three states |0) = \Nl,Nr), \L) = \Nl + 1,Nr), 
and \R) = \Nl,Nr + 1), which correspond to many- 
particle ground states with Nl{+1) electrons in the left 
and Nr{+1) electrons in the right dot. The correspond- 
ing ground state energies El of \L) and en of \R) are 
in the window /i^ > £lt£r > PR- We assume that the 
Coulomb charging energy Uc is the largest energy scale, 
and it is not possible to charge the double dot with more 
than one additional electron. In Q, no enhanced tunnel 
current was observed for e := — e^^ < at low tem- 
peratures so that excited many-body states play no role. 
In particular, Uc ^ 1 meV was one order of magnitude 
larger than the external source drain voltage Vsd- This 
situation has to be contrasted with the regime Vsd ^ Uc 
(absence of blockade effects, fl^). 

We define operators ul = \L){L\, ur = \R){R\, p = 
\L){R\, SL = |0)(L|, sr = \0)(R\, and the total system 
Hamiltonian H as the sum of the dot, the phonon, the 
reservoir and the electron-phonon interaction 
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Here, the tunneling between left and right dots is de- 
scribed by a single tunnel matrix element Tc- In the 
standard tunnel Hamiltonian Hy , Vk and Wk couple the 
dot to a continuum of channels k of the left and right 
electron reservoir Hres- The spin of the electron does 
not play any role here and is suppressed. The term Hp 
describes the lattice vibrations in harmonic approxima- 
tion; the creation operator for a phonon of mode Q is 
aq . The electron-phonon matrix elements are defined by 
aq AQ(L|e'Q"-|L) and /3q := XQ{R\e'^''\R) , where Aq 
is the matrix element for the interaction of 2DEG elec- 
trons and phonons. We have already neglected the non- 
diagonal term of the electron- phonon interaction which 
contains matrix elements 7q := AQ(_L|e**^''|i?). We have 
checked in a seperate master equation calculation that 
such terms modify only weakly the tunnel current and 
do not lead to the oscillatory phenomena observed in Q , 
which is due to the non-perturbative shake-up process 
that we describe in the following. 

Suppose an electron tunnels between two regions of 
space (L and R) and interacts with a phonon field. With 
the interaction of the form Ha/s, Eq.(|l|), the electron- 
phonon coupling locally changes the energy of the elec- 
tron, depending on whether it is in L or in R. During 
tunneling, its wave function experiences an additional a 
phase shift e** due to this coupling. This phase shift is 
zero if the coupling is identical in both regions, i.e. if 
Q!Q = /3q. If, as we show below, aq ^ /3q, the phase 
shift depends on the state of the phonon field. From 
the point of view of the phonons, initial (before the tun- 
neling) and final (after tunneling) phonon states are no 
longer the same. From the point of view of the electron, 
the effective tunnel amplitude changes non-trivially. In 
particular, it becomes time-dependent. 

We introduce now a unitary polaron transformation 
psf of the Hamiltonian that naturally leads to the phase 
factors mentioned above. For any operator O, we de- 
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ulA + ub.B with A 



Eq(1/^q) ("Qflq - a-Qflq) and B := Eq(l/^Q) 
(/SqAq — /^-qflq). This leads to renormalized energies 

el ^ eL~ Eq I"qIV'^Q and e]^^ er- J2q I/3qIV^Q> 
and a renormalized intra-dot tunneling Hamiltonian 
Ht — Tc{pX -I- p'I'X'l'). Here, the phase operator X — 
riq ((o^Q ~ /5q)/i^q) is the product of unitary dis- 
placement operators Dq{z) := exp(zaQ — 2*aq), where 
the operation of Dq{z) on the boson vacuum creates a 
coherent state of the boson field mode Q. The factors X 
and X'^ in the tunnel Hamiltonian Ht drastically change 
the transport properties of the double dot. 

We assume that the coupling to the left and right elec- 
tron reservoirs is weak such that a standard Born and 
Markov approximation holds. In contrast to this, we 
calculate to all orders of the miradot tunneling Tc and 



the electron-phonon coupling because the renormaliza- 
tion of the tunneling by the phase factors X is a non- 
perturbative effect. From the Liouville equation for the 
total density matrix of the system, one obtains the equa- 
tions of motion (cp. pl] , p^ ) 
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where F^ := 27r j:^V^ 5(sl - e^), F^ := 2^ W^k 
S{eR — e^), and the chemical potential (pr) of the 
left (right) electron reservoir is assumed to be far above 
(below) El (er) so that no electrons can tunnel from the 
left dot into the left reservoir, or from the right reservoir 
into the right dot. Furthermore, e = £l — Er, where 
we neglect the difference in the energy renormalization 
in both dots, p{t) — pe"^*Xt, and Xt denotes the time- 
evolution of X with Hp. 

We assume the phonon system to be in thermal equi- 
librium at all times. We use a decoupling of the reduced 
density matrix p{t') of the dot, p{t') « Pph'^'^ph{p{^')}- 
We introduce the Laplace transform 
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and nL{z) dte~^^ {nL)t etc. for z > 0. By trans- 

forming Eq.(|^) into z-space, the expectation value of the 
current operator / := iTc{p — p^) (we set the electron 
charge e — I for convenience), can be obtained in the 
stationary limit t ^ oo from the l/z-coefficient of the 
/(z)-expansion into a Laurent series for z ^ 0. The re- 
sult is 
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Here, — C^{5 0). One can verify that in the limit 
of no phonon coupling, = i/e, thus obtaining a sim- 
ple Lorentzian curve for the stationary tunnel current 
as a function of e [O. The appearance of Tc in the 
denominator of Eq.(H) indicates that this result is non- 
perturbative, i.e. valid to all orders in Tc- The modifi- 
cation of this curve by the electron-phonon interaction 
is completely described by the form of which we will 
discuss now. 

For a harmonic phonon system in thermal equilib- 
rium at temperature T, one finds the correlation function 
{XtXl)o = e-*(*-*') with 
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The matrix elements aq and /3q are aq = 
Aq/ d3xe'Q''pL(x), /3q = Aq/ rf^xe^Q'^pflCx), with the 
local electron densities pl(x) and Pfl(x) in the left and 
the right dot, respectively. The exact form of both pl(x) 
and /9_r(x) depends on the shape of the dots and on 
the number of electrons and Nn. In the stationary 
state for f — > cx), both densities can be assumed to be 
smooth functions of x, ~ pe(x — x^), P_r(x) ~ 

Pe(x — Xfl), where we have assumed that both left and 
right electron densities are described by the same profile. 
Here, Pe(x) is relatively sharply peaked around zero. One 
obtains aq = AQe*'^'-^P(Q), /3q = AQe'Q'-«P(Q), where 
P(Q) = / d^xe^Q^^pelx). Therefore, /3q = oqc'^'^ with 
d = rii — TL. The matrix elements aq and /3q thus coin- 
cide up to the phase factor e**^'*. This phase factor is the 
origin of the oscillations in the effective density of states 
p{uj) of phonon modes that couple to the tunnel process. 
We note that the explicit form of aq and /3q can be cal- 
culated exactly when the two states \L) and \R) are few 
particle wave functions. 

For sharp charge density profiles with Fourier trans- 
form P(Q) ^ 1, we find 
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where ujd '■— Cs/d, d — |d|, Cs is the longitudinal speed 
of sound, and a smooth exponential cutoff lOc is an ef- 
fective Debye frequency. The latter is due to the finite 
extension / of the wave functions in the dots leading to a 
cutoff of phonons with wave vectors Q = |Q| > l/l. The 
phonons are assumed to be piezoelectric acoustical bulk 
(three-dimensional) modes with an interaction matrix el- 
ement |AqP = /VcsQ, A^ = hP/2pM, where wq — CgQ 
is the phonon dispersion, V the volume and pM the mass 
density of the crystal, and P the piezoelectric coupling 
constant. Dynamical screening effects of the 2DEG have 
been absorbed into the value of P. Using typical GaAs 
parameters we obtain g « 0.05. 

The most important feature of p{uj), Eq.(||), are the 
oscillations on the scale LUd = Cg/d, (Fig. ||, inset). In 
fact, these lead to the oscillations in the current pro- 
file due to phonon emission in the experiment Q: Fig. |l| 
shows the stationary current Eq. (^) at different tempera- 
tures for parameters close to the ones of the experiment. 
With d = 200 • 10-^m and = 5000to/s, we obtain 
fiud = IQ-bp-eV, which is the scale on which the oscil- 
lations occur in Q. The cutoff energies are assumed as 
hujc — 1 meV and hd — 1/icV. A small but finite value of 6 



in Cg {S O) simulates off-diagonal electron-phonon pro- 
cesses which are not included in p{uj) for hu) = e 0. At 
low temperatures, there is a broad oscillatory shoulder for 
e > (spontaneous phonon emission). It reflects the os- 
cillations in p{uj) which in turn determines Cg in Eq. (|^). 
Its real part 5RCe is proportional to the probability den- 
sity for inelastic tunneling from the left to the right dot 
with energy transfer e ||^. In the limit = 0, one finds 
= {2Tr/T{g))eS-^e-''9{e) at T = 0, where only spon- 
taneous phonon emission is possible. This shows that the 
effect is non-perturbative in the electron-phonon coupling 
9- 
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FIG. 1. Stationary tunnel current, Eq.(^), as a function of 
the energy difference e between left and right dot ground state 
energies. Dimensionless electron-phonon coupling parameter: 
g = 0.05. Inset: effective density of states p{u>) of phonon 
modes, Eq.(|) and Eq.(|). 

As an important consequence, even for small coupling 
constants g and T ^ 0, the electron-phonon interaction 
can dominate the non- linear electron transport. Further- 
more, we find an energy dependence of the current for 
e ^ 50 /xeV between 1/e (larger g) and (smaller 
g), depending on g. Note that the tunnel current / is 
determined by the complex function which contains 
all effects of energy renormalization due to the electron- 
phonon coupling. At higher temperatures, the current 
at the absorption side e < increases faster than at the 
emission side where the oscillations start to be smeared 
out. For e < a shoulder-like structure appears, consis- 
tent with the observation in Q. 

We note that our results are consistent with a previous 
model for inelastic tunneling via pairs of impurities. 
There, perturbation theory in the electron-phonon de- 
formation potential coupling was used to obtain the non- 
linear current (without oscillations) for large area tunnel 
junctions after averaging over an ensemble of impurity 
pairs. 
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The spontaneous emission of phonons in double dots is 
closely analogous to the subradiant spontaneous decay by 
emissions of photons from a laser-trapped two-ion system 
1^. There, the interaction of the atomic dipoles d^, at 
positions r.^, i — 1,2 with the electromagnetic field has 
the form ^ 

Heph = ^gQs(diexpiQri +d2CxpiQr2), (7) 

Qs 

with gQs = —i{2TTcQ/Vy^^eQ,s with the speed of light c, 
and polarization vectors eq.s for polarization direction s. 
This leads to rates for super- and subradiant decay of the 
form 11 r(g)± = roiQ)[l±asm{Qd)/{Qd)] with a = 1 
(a = 3/2) if the vector character of the light is (not) 
neglected, Q = lu/c, and ro{Q) oc Q^. The subradiant 
channel (— ) is due to the decay of the singlet state that 
corresponds to the difference di exp iQrj^ — d2 exp iQr2 
of the dipole moments. The interaction with a phonon of 
mode Q in the double dot is oc exp iQ,Vi+nR expiQrj. 
Thus, it has the same 'interference form' as in the two- 
ion case, Eq.(Q). The tunnel current is modified by the 
phase difference of the electron before and after tunneling 
whence piuj), Eq.(|^), corresponds to the subradiant rate 
r_ . Although the microscopic mechanism is different in 
both cases (for light the rates Fq are oc w"^, for piezoelec- 
tric phonons oc the interference term sm{Qd) / Qd 
in both cases is due to the 'interference of matrix ele- 
ments' in the interaction Hamiltonian and has the same 
physical origin. 

In the atom-trap experiment, the Dicke effect, i.e. the 
existence of two different radiation channels T±, has been 
verified by changing the ion distance. The experimen- 
tal data of the double quantum dot seem to indi- 
cate that for larger distance d of the dots the oscillations 
oc sai{e/fiujd) as a function of e become faster which is 
consistent with ujd = Cs/d in Eq.(^. 

In conclusion, we have found an interference effect 
which explains the recently observed phonon emission 
spectrum in the transport through coupled quantum 
dots. The analogy between artificial and real multiple 
atom systems allows to speculate about further realiza- 
tions of coherent optical effects in semiconductor quan- 
tum dot experiments, similar to the oscillatory Dicke ef- 
fect we have found recently jl^. In the case of phonons 
discussed here, the coherence (oscillations due to interfer- 
ence) showed up within a dissipative process, i.e. spon- 
taneous emission. In order to further study such 'coher- 
ent dissipative phenomena', we suggest systems of two 
of more double quantum dots which are coupled via the 
common phonon field. The non-linear transport proper- 
ties are then determined by the electron-phonon coupling 
both within each and between the double dots. 
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